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Method of determining impulse responses from a medium 
in relation to the transmission of waves between 
different points 

5 The present invention concerns methods of determining 
impulse responses from a medium in relation to the 
transmission of waves between different points. 



More particularly, the invention relates to a method of 
10 determining impulse responses of a medium in relation 
to the transmission of waves between different points, 
method comprising: 

(a) at least one step of emission in the course of 
which waves are emitted into the medium by generating 

15 signals ei(t) on the basis of a number N of emission 
points belonging to the medium, where N is an integer 
at least equal to 2 and i is an index lying between 1 
and N which designates one of said N emission points, 

(b) at least one step of reception in the course of 
20 which signals rj (t) are picked up from said waves after 

transmission in said medium, at a number M of reception 
points belonging to the medium (these reception points 
may be coincident with the emission points) , where M is 
a non-zero natural integer and j is an index lying 
25 between 1 and M which designates one of said M 
reception points, 

(c) and at least one step of determination of said 
impulse responses hij (t) between each emission point i 
and each reception point j on the basis of the signals 

30 emitted ei(t) and picked up rj (t) . 



The impulse responses may be determined for the 
purposes of focusing waves in the medium, for example 
for the purposes of imaging the medium or of selective 
35 communication between various points of the medium, in 
particular within the framework of the implementation 
of the "D.O.R.T." procedure (Decomposition de 
l'Operateur Retournement Temporel [Decomposition of the 
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Time Reversal Operator] ) described in particular by 
Prada et al., (C. Prada and M. Fink, "Eigenmodes of the 
time reversal operator: a solution to selective 
focusing in multiple-target media." Wave Motion, 20, 
5 pp 151-163 (1994) ) . 

The process of acquisition of the impulse responses 
between the points in question of the medium is, 
however, relatively slow, insofar as it is performed 
10 sequentially, for the various emission points i one 
after the other, as described for example in the 
document WO-A- 02/32 316 within the framework of a 
method of focusing of acoustic waves other than the 
abovementioned "D.O.R.T." procedure. 

15 

This slowness lengthens the time of implementation of 
the method, and when the medium is changeable, it does 
not make it possible to obtain the impulse responses 
quickly enough to be able to be used validly thereafter 

2 0 in this medium. 

The present invention is aimed in particular at 
alleviating these drawbacks. 

25 For this purpose, according to the invention, a method 
of the kind in question is characterized in that in the 
course of step (a) , said N emission points are made to 
simultaneously emit the signals ei(t), these signals 
ei(t) having a duration T and each being a sum of n 

3 0 substantially monochromatic elementary signals, of like 

amplitude and of respective frequencies f 0 ,i+k.Sf, where 
f 0 ,i is a predetermined eigenf requency at the point i, k 
is an integer lying between 0 and n, n is an integer at 
least equal to 2 and Sf is a predetermined frequency 
35 interval, the respective eigenf requencies f 0 ,i at the 
various points i being distinct from one another and 
lying in a frequency band of width 8f , 
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and in that in the course of step (c) , each impulse 
response hij (t) is calculated on the basis of a signal 
of correlation between the signal ei(t) emitted at the 
point i and the signal rj (t) picked up at the point j . 

5 

By virtue of these provisions, the impulse responses 
between N emission points and M reception points 
(possibly coinciding with the emission points) are 
calculated simultaneously, thereby accelerating the 
10 process for acquiring the impulse responses. 

In various modes of execution of the method according 
to the invention, recourse may possibly be had 
furthermore to one and/or other of the following 
15 provisions: 

the respective eigenf requencies f 0/ i at the various 
points i are separated pairwise by an offset 8f/N; 

in the course of step (c) , said correlation signal 
is windowed by means of a gate function 7t(t) of width 
20 l/8f; 

in the course of step (c) , the impulse responses 
hij (t) are determined through the formula: 
hij(t) = U(t) .Jei(0 - t) .rj(0)d0; 

the waves transmitted in the medium between the 
25 emission points and the reception points are acoustic 
waves; 

in the course of step (a) , the medium where the 
waves are emitted is reverberant; 

the frequency interval Sf is less than or equal to 
3 0 1/t, where t is the temporal dispersion of the medium 
(that is to say the duration of the signal picked up 
after emission of an impulse signal into the medium) ; 

the frequency interval 5f is substantially equal 
to 1/t, where x is the temporal dispersion of the 
3 5 medium; 

the duration T is at least equal to N/Sf ; 

the duration T is at least equal to N.x, where t is 
the temporal dispersion of the medium; 



WO 2004/086557 - 4 - PCT/FR2004/000684 

" the elementary signals exhibit random phases; 

the waves are emitted with a certain passband, the 
frequencies fOi comprise a minimum frequency fO and the 
number n is determined so that the frequency band lying 
5 between fO and f 0+ [ (n+1) .8f ] substantially overlaps 
said passband; 

the reception points are coincident with the 
emission points. 

10 Other characteristics and advantages of the invention 
will become apparent in the course of the following 
description of two of its modes of execution, given by 
way of non- limiting example, with regard to the 
appended drawings . 

15 

In the drawing, figure 1 is a diagrammatic view of a 
device suitable for implementing the present invention, 
in a mode of execution. 

2 0 In the non- limiting example represented in the drawing, 

the method according to a mode of execution of the 
invention is implemented by a device 1 which is 
suitable for emitting and receiving waves in a medium 
2, homogeneous or otherwise, generally reverberant, 
25 which may be solid, liquid or gaseous. The waves in 
question are preferably acoustic and the description 
hereinbelow will refer only to such acoustic waves, but 
said waves could possibly be of some other nature, for 
example of electromagnetic nature. 

30 

The device 1 comprises for example a computer 3, or any 
other similar control device, which controls the 
emissions of acoustic waves and the listening thereto 
in the medium 2. For this purpose, the computer 3 may 

3 5 for example control a dedicated electronic device 4 

which may in particular comprise an electronic central 
unit CPU associated with a central memory MEM, this 
central unit CPU itself controlling N pathways (N being 
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an integer at least equal to 2) each linked to a 
transducer Tl-TN making it possible to emit and to pick 
up acoustic signals in the medium 2. All the pathways 
are synchronized with one another. 

5 

Each measurement pathway of the device 4 may comprise 
for example a sampler El -EN linked to the corresponding 
transducer Tl-TN, and each sampler may itself be linked 
to a memory associated with a central unit, referenced 
10 jointly Cl-CN. 

The device 1 may be used in particular for the purposes 
of imaging the medium 2, or else to establish a 
selective communication with a point of said medium, or 

15 the like, for example using the "D.O.R.T." procedure 
(Decomposition de l'Operateur Retournement Temporel) 
described in particular by Prada et al . , (C. Prada and 
M. Fink, "Eigenmodes of the time reversal operator: a 
solution to selective focusing in multiple- target 

20 media." Wave Motion, 20, pp 151-163 (1994)). 

In all cases, it is necessary to determine the impulse 
responses between the transducers Tl-TN. 

25 For this purpose, the computer 3 controls the device 1 
in such a way that it follows successively the 
following steps: 

(a) at least one step of emission in the course of 
which waves are emitted into the medium by generating 

30 signals ei(t) on the basis of the N transducers Tl-TN 
(these N transducers may as appropriate constitute only 
part of the transducers linked to the device 4, in 
which case the impulse responses pertaining to the 
whole set of transducers would be acquired in several 

35 successive passes, these passes nevertheless remaining 
much fewer in number than the total number of 
transducers) , i being an index lying between 1 and N 
which designates one of said N transducers; 
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(b) at least one step of reception in the course of 
which the N transducers Tl-TN pick up signals rj (t) 
from said waves after transmission in said medium, 
where j is an index lying between 1 and N designating 

5 one of said transducers; 

(c) and at least one step of determination of said 
impulse responses hij (t) between each emission point i 
and each reception point j on the basis of the signals 
emitted ei(t) and picked up rj (t) . 

10 

To make it possible to acquire in one go all the 
signals necessary for the determination of the impulse 
responses relating to the N abovementioned transducers 
(that is to say all or part of the transducers linked 

15 to the device 4, as explained hereinabove), all the 
signals ei(t) are emitted simultaneously in the course 
of step (a) by the N transducers, but these emitted 
signals are orthogonal to one another so that the 
information specific to each emission point may 

20 thereafter be separated from the other information in 
the signals picked up rj (t) . 

Thus, each emitted signal ei(t) is equal to the sum of 
n substantially monochromatic elementary signals, of 
25 like amplitude and of respective frequencies f 0 ,i+k.8f, 
where : 

f 0 ,i is a predetermined eigenf requency specific to 
transducer i, 

k is an integer lying between 0 and n, 
30 n is an integer at least equal to 2, 

and 5f is a predetermined frequency interval, 
advantageously less than or equal to 1/x and preferably 
substantially equal to l/x, where x is a mean value of 
the temporal dispersion of the medium 2, that is to say 
35 the duration of the signal picked up after emission of 
a pulsed signal (Dirac function) in the medium 2 . 
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The respective eigenf requencies f 0 ,i specific to the 
various points i are distinct from one another and lie 
in a frequency band of width 5f : advantageously these 
frequencies f 0 ,i are separated pairwise by an offset 
5 8f/N. 

Moreover, if BP is the passband with which the acoustic 
waves are emitted and received by the transducers, and 
if the lowest of the frequencies fOi is dubbed fO, then 
10 fO and n may advantageously be determined so that the 
band of the frequencies of the elementary signals (band 
lying between fO and f 0+ [ (n+1) .8f ] ) substantially 
overlaps said passband. 

15 Advantageously, the elementary signals making up the 
signals ei(t) exhibit random phases. 

All the signals ei(t) have the same duration T, which 
is advantageously at least equal to N/5f, and 
20 preferably at least equal to N.t. 

In the course of step (c) , each impulse response hij (t) 
is calculated on the basis of a signal of correlation 
between the signal ei(t) emitted at the point i and the 
25 signal rj (t) picked up at the point j. 

In the example considered here, this correlation signal 
equals : 

cijit) = jei(0) .rj(t - 0)dO (this integral is done for 

30 example over a duration at least equal to T, the start 
of which coincides with the start of the signal rj (t) ) . 

Advantageously, this correlation signal is windowed by 
means of a gate function 7i(t) of width l/5f (which 
35 equals 1 during a period of l/5f and 0 outside of this 
period) , in which case the impulse responses hij (t) may 
advantageously be determined through the formula: 
hi jit) = Yl(t) .cij(t) . 
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It will be noted that, in the particular example 
represented in the drawing and described hereinabove, 
the reception points j at which the acoustic wave is 
5 picked up are coincident with the emission points and 
constituted by the transducers Tl-TN, but the reception 
points in question may as appropriate be constituted by 
a second array of transducers (not represented) , 
distinct from the array of transducers Tl-TN. 

10 

In the general case, the acoustic waves which are 
generated by the emission of the signals ei(t) are 
therefore picked up at M reception points, M being a 
non-zero integer, which M reception points may either 
15 be coincident with the N emission points (in which case 
M = N) , or not be completely coincident with the 
emission points, or be totally distinct from the 
emission points. 



